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Abstract 

Signals consisting of a sequence of pulses show that inherent origin of the 1/f 
noise is a Brownian fluctuation of the average interevent time between subsequent 
pulses of the pulse sequence. In this paper we generalize the model of interevent 
time to reproduce a variety of self-affine time series exhibiting power spectral density 
S(f) scaling as a power of the frequency /. Furthermore, we analyze the relation 
between the power-law correlations and the origin of the power-law probability dis- 
tribution of the signal intensity. We introduce a stochastic multiplicative model for 
the time intervals between point events and analyze the statistical properties of the 
signal analytically and numerically. Such model system exhibits power-law spectral 
density S(f) ~ 1/f 13 for various values of /3, including /3 = 1/2, 1 and 3/2. Explicit 
expressions for the power spectra in the low frequency limit and for the distribution 
density of the interevent time are obtained. The counting statistics of the events is 
analyzed analytically and numerically, as well. The specific interest of our analysis is 
related with the financial markets, where long-range correlations of price fluctuations 
largely depend on the number of transactions. We analyze the spectral density and 
counting statistics of the number of transactions. The model reproduces spectral 
properties of the real markets and explains the mechanism of power-law distribution 
of trading activity. The study provides evidence that the statistical properties of the 
financial markets are enclosed in the statistics of the time interval between trades. 
A multiplicative point process serves as a consistent model generating this statistics. 

Keywords: Stochastic processes; Econophysics; Financial markets; Point pro- 
cesses; 1/f noise 
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1 Introduction 

Complex collective phenomena usually are responsible for the power-laws which are 
universal and independent of the microscopic details of the phenomenon. Examples in 
physics are numerous. Power-laws are intrinsic features of the economic and financial 
data, as well. The aim of this contribution is to analyze a relationship between the 
origin of the power-law distribution and the power-law correlations in financial time 
series. 

There are numerous studies of power-law probability distributions in various 
economic systems [1-15]. The key result in recent findings is that the cumulative 
distributions of returns and trading activity can be well described by a power-law 
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asymptotic behavior, characterized by an exponent A ~ 3, well outside the stable 
Levy regime < A < 2 [13-15]. Empirical studies confirm the power-law asymp- 
totic behavior with stable exponent for different time scales. Since these empirical 
distributions are neither Levy stable nor invariant under addition, new stochastic 
models with long range correlations and power-law asymptotic behavior are of great 
interest [14]. 

The random multiplicative process built into the model of wealth distribution 
yields Pareto power-law [2]. The generalized Lotka-Volterra dynamics is in the use 
for various systems including financial markets [9]. However, these models generically 
lead to non-universal exponents and do not explain the power-law autocorrelations 
and observable spectral density of financial time series [11]. The time-correlations in 
the financial time series are studied extensively as well [13-17]. Recent investigations 
[13-16] provide empirical evidence that the long-range correlations for volatility are 
due to the trading activity, measured by a number of transactions N. Therefore, 
various stochastic models of trading activity or for waiting time between transactions 
have to be analysed for reproduction the statistics of volatility in financial markets. 

Recently we tried to adapt the model of 1/f noise based on the Brownian mo- 
tion of time interval between subsequent pulses, proposed in [18-21], to model the 
share volume traded in the financial markets [22]. The idea to transfer long time 
correlations into the stochastic process of the time interval between trades or time 
series of trading activity is in consistence with the detailed studies of the empirical 
financial data [13, 16] and can reproduce the spectral properties of the financial time 
series [22]. However, the investigation of the model revealed, that the simple additive 
Brownian model of time interval between trades failed to reproduce the power-law 
probability density function (pdf) of the trading activity [23]. 

On the other hand, several authors have shown empirically that the fluctuations of 
various financial time series possess multifractal statistics [24-28]. Therefore, in this 
paper we introduce the stochastic multiplicative model for the time interval between 
trades and analyze the statistical properties of the trading activity analytically and 
numerically. The model is in consistence with the results of statistical analysis of 
the empirical financial time series and reveals the relationship between the power-law 
probability distribution and the power-law spectral density of the financial time series. 
Moreover, the analytically solvable model [18-21] results only in exact 1/f noise while 
in the papers [22, 23] essentially only numerical analysis has been undertaken. Here 
we generalize the model defining the time series exhibiting the power spectral density 
1/f 13 with 0.5 < (3 < 1.5 and the power-law probability density of the trading activity. 
Explicit analytical expressions for the power spectra and for pdf of the trading activity 
N are obtained and analysed. 

2 The model 

We consider a signal I(t) as a sequence of the random correlated pulses 

I{t) = ^2a k 5(t-t k ) (1) 

k 

where a k is a contribution to the signal of one pulse at the time moment t k , e.g., a 
contribution of one transaction to the financial data. Signal represents a point 
process used in a large variety of systems with the flow of point objects or subsequent 
actions. When a k = a is constant, the point process is completely described by 
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the set of times of the events {t^} or equivalently by the set of interevent intervals 
{jk = ^fc+i ~~ ifc}- Various stochastic models of can be introduced to define a 
stochastic point process. In papers [18-21] it has been shown analytically that the 
relatively slow Brownian fluctuations of the interevent time yield 1/f fluctuations 
of the signal 

Power spectral density of the signal I(t) is 



S(f) 



lim / — — 

T^oo\ T 



2a 



I(t) exp(-i2irft)dt 



(2) 



2&2o\"r" ^ S exp(-i2vr/A(A : ;g)) 



fc ^min 9 ^min ^ 



where the brackets (...) denote the averaging over the realizations of the process, 
T = tf — ti is the observation time, A(k; q) = tk+ q — t k is the time difference between 
the pulses occurrence times tk +q and t^, while k m i n and /c max are minimal and maximal 
values of the index k in the interval of observation T. 

We adapt and generalize the model of 1/f noise previously proposed in [18-21] 
to model the share volume traded in the financial markets [13, 16]. It is useful to 
define a discrete time series with equal time intervals Td- Integration of the signal 
I(t) over subsequent time intervals of length results in a discrete time series and, 
by analogy with financial time series, we will call it the volume Vj, 



tj+r d 



I(t)dt 



^2 

tj<t k <tj+T d 



JTd- 



(3) 



The number of trades Nj in the time interval [tj,tj + Td] is defined by the same 
equation © with a& = 1 and Nj = Vj. Then the power spectral density S(f s ) of the 
discrete signal Nj may be calculated by the FFT as 



S(fs) = — 
\ T d n 



^iVxpH27r( S -l)(i-l)) 

3=1 



(4) 



Here the discrete frequencies f s = (s — 1)/T, s = 1, 2, ...n and T = Tdn. For frequen- 
cies f s <C 1/Td the power spectral density defined by equation @ coincides with that 
given by equation p|. 

In this paper we investigate the statistical properties of the time series Nj, when 
the sequence of interevent time Tk is generated by a multiplicative stochastic process. 
First of all the multiplicativity is an essential feature of processes in economics [1,7]. 
Multiplicative stochastic processes yield multifractal intermittency and are able to 
produce power-law probability distribution functions. We base our study on the 
generic multiplicative process for the interevent times, written as 



Tk+l =T k + lT k 



2/j,- 



1 . a 



(5) 



Here the interevent time T k fluctuates due to the external random perturbation by 
a sequence of uncorrelated normally distributed random variable with zero ex- 
pectation and unit variance, a denotes the standard deviation of the white noise and 
7 <C 1 is a damping constant. 
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We will restrict the diffusion of the interevent time according to equation (JSJ to 
the finite interval [T m i n ,r max ], i.e., < r m ; n < < r max . The most simple case 
is the pure multiplicativity of i.e., when fj, = 1. Other values of /x reproduce 
the power-laws as well and explicit expressions can be derived without the loss of 
generality. 

The iterative relation © can be rewritten as Langevine stochastic differential 
equation in fc-space. 

^=^- 1 +r^(k). (6) 

Here we interpret k as continuous variable while 

(H(k)H(k')) = 5(k-k'). 

The steady state solution of the corresponding stationary Fokker-Planck equation 
with zero flow gives the probability density function for r k in the /c-space (see, e.g., 
[29]) 

Pk(r k ) = Ct£, a = 2 7 /cx 2 - 2/x (7) 
and C has to be defined from the normalization 

/'''"max 

/ P k (T)dr = l. 

J 7"min 

The solution Q assumes Ito convention involved in the relation between expres- 
sions ©, © and (J7J). 

As the probability distribution function Q follows a power-law, we expect that 
this would result in the power-law behavior in other statistics, as well. In the next 
Section we will obtain the power spectral density defined by equation (|2|). 



3 Power spectral density 

The power spectral density is a well-established measure of long time correlations 
and is widely used in stochastic systems. As it has been already shown [18-21] the 
point process with the Brownian interevent time exhibits 1// noise. 

We will derive the formula for the power spectral density of the multiplicative 
stochastic point process model, defined by equations © and © for the interevent 
time. 

Let us rewrite equation (J2J) as [21] 

In 2 

5(/)= T lim o ^^XA( fc:g )(2vr/) (8) 

k,q 

where XA(fc:g)(27r/) is the characteristic function for the probability distribution of 
A(k;q). For the normal distribution of A(k;q) formula ((HJ) takes the form 

2tj2 

S(f) = lim — V exp{i2vr/ <A(fc; q)) - 2^fo\{k- q)} (9) 

1 — >oo 1 * — • 

k,q 

where cr^ is the variance of A(fc; q), the time difference between the pulses occurrence 
times tk+q and t k - For /i = 1 the time difference A(k; q) may be expressed from the 
solution of multiplicative stochastic equation (jHJ) as 

q I 

A(k; q) = Y, T k exp{( 7 - -a 2 )l + a £ e 3 }. (10) 
i=i j=i 
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Averaging over the normal distributions of uncorrelated e~ results in explicit ex- 
pressions for the mean (A(k;q)) and variance a\(k;q). 

In general, for any //, the perturbative solution of equation (jHl) yields 

(A(k;q)) = r.q + lr^q 2 + o( 7 2 ), (11) 



<r\{k-q) = rjj'*- 



3 



In the low frequency limit / <C t^ 1 we can replace the summation over k and q 
by the integration and take into account only the first order terms of equation (jlljl 
in the expression for the power spectral density @. This yields 

S,(f) = ^Re J drr a J 'exp{z27r/(r g + ^-Y)}d g , 



' mm 



2Ca 2 / 7 \ 3 - 2 m 

^erf (■^minj 2-maxj) (1^) 



^raax 

4rf(^min,^max) = Re J exp j-i - ^ J | erfc(V -ix)x^^~^ dx (13) 



where 

^ = in) 



h — k ■ 
Tnax "rnin 



is the expectation of r^. Here we introduce the scaled variable and 



T min > X max - 



Expression of the power spectral density (|12j) is appropriate for the numerical 
calculations of the generalized multiplicative point process defined by equations Q 
and ©. In the limit r m i n — ► and T max — > oo we obtain an explicit expression for 



V / 7TT(3 - 2y«)/ \TTfJ cos ( 2(3-^) ^ 

Equation (|14() reveals that the multiplicative point process (jSJ results in the power 
spectral density £(/) ~ l//' 3 with the scaling exponent 

Let us compare our analytical results (|12)1. (|13j) and (|14j) with the numerical 
calculations of the power spectral density according to equations © and (JSJ). In 
Fig. [l\we present the numerically calculated power spectral density S(f) of the signal 
I(t) for fi = 1 and a = 27/a 2 — 2 = 0, -0.4 and +0.5. Numerical results confirm 
that the multiplicative point process exhibits the power spectral density scaled as 
S(f) ~ 1//^- Equation (|12|) describes the model power spectral density very well in 
a wide range of parameters. The explicit formula (|14|) gives a good approximation of 
power spectral density for the parameters when ~ 1. In Fig.Elwe present numerical 
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Figure 1: The power spectral density S(f) of the signal I(t) versus frequency / calculated 
from the model described by equations (JTJ), (j2J and (JHJ) with jjl = 1. The sinous curves 
represent the results of the numerical simulations averaged over 5 realizations of 100000 
impulse sequences with a = 0.1, the straight lines represent the analytical approximation 
given by integral formula (fl~2|) with r m i n = 10 -6 ; r max = 1 and the dashed lines show results 
of explicit formula JH). (a) 7 = 0.008, (b) 7 = 0.01 and (c) 7 = 0.0125. 

results of the power spectral density calculated for the multiplicative point process 
with n = 0.5. In this case analytical expressions (|12j) . (|13|) and 1)14(1 describe power 
spectral density very well in a wide range of the parameter 7. These results confirm 
the earlier finding [18-21] that the power spectral density is related to the probability 
distribution of the interevent time and 1/f noise occurs when this distribution is 
flat, i.e., when a = 0. 

It is likely that such a stochastic model with parameters in the region 0.5 < (3 < 
1.5 may be adaptable for a wide variety of different systems. In this paper we will 
investigate applicability of the model for the financial market. 

4 Power-law probability distribution function 

The statistics of the trading activity N, i.e., the number of trades per time interval 
Trf is the major task of this paper. The consistency of any model of financial markets 
depends on its ability to reproduce empirically measured probability distribution of 
the number of trades N, defined by equation (j3J). Intuitively, the pdf of ./V is related 
with pdf of the rate v = 1/r and with pdf of r. The pdf of in A;-space is given by 
expression (JTJ). It is obvious that in the actual time t the pdf of r may be written as 

P t (r) = P k (r)r/f = C'r a+1 (16) 
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Figure 2: The same as in Fig. □ but for // = 0.5. (a) 7 = 0.001, (b) 7 = 0.0025, (c) 
7 = 0.005 and (d) 7 = 0.0075. 

where C = C/f is a new normalization constant. 

For the pure multiplicative model with \x = 1 equations (|10|) and (jllj) define a 
relationship between iV and r after the substitution (A(k;q)) — > r^, — > r and 
Q ^ N, i.e., 

r d = r7V+^r7V 2 . (17) 

This relationship may be used for definition of pdf P(N) from the relation 
P{N)dN = Pt{r)dT. We have 

CV d 2+a (l + 7 iV) 
P(iV) ~ AT3+a(i + 2iV)3+a- ( 18 ) 

This equation yields limiting cumulative distributions for N, 

+00 

P > (N)= [p(N)dN~{ (19) 

AT 

Let us compare probability distribution function for ./V obtained from the numer- 
ical simulation of the model with the analytical results ()18|) and 1)19|) . It is convenient 
to use the amplitude versus rank plot (see, e.g., [30] ) for definition of the exponent 
A of the cumulative distribution P> (N) ~ 1 /iV A . We plot the value of the n-th rank 
variable N(n) as a function of its rank n. If P > (N) behaves asymptotically when 
n — > as a power-law, one gets a straight line in log-log coordinates of N and n, 
with the slope equal to — 1/A. In Fig. El we present the plot for N(n) calculated 
numerically for pure multiplicative model (/u = 1) with the parameters yielding the 
best fit to the empirical data. 
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Figure 3: Amplitude versus rank, N(n), plot defining cumulative distribution P>(N) ~ 
1/N X with the parameters of the model \x = 1, a = 0.1 and 7 = 0.0085. Dashed line 
represents asymptotic plot N(n) ~ n _1 / A with A = 3.4. 

5 Discussion and conclusions 

We have introduced a multiplicative stochastic model for the time intervals between 
events of point process. Such a model of time series has only a few parameters 
defining the statistical properties of the system, i.e., the power-law behavior of the 
distribution function and the scaled power spectral density of the signal. The ability 
of the model to simulate 1// noise as well as to reproduce signals with the values 
of power spectral density slope (3 between 0.5 and 1.5 promises a wide variety of 
applications of the model. 

Let us present shortly the possible interpretations of the empirical data of the 
trading activity in the financial markets. With a very natural assumption of transac- 
tions in the financial markets as point events we can model the number of transactions 
Nj in equal time intervals as the outcome of the described multiplicative point 
process. We already know from available studies [13] that the empirical data exhibit 
power spectral density in the low frequency limit with the slope j3 ~ 0.7. For the pure 
multiplicative model fi = 1 this corresponds to the case when a = 2j/a 2 — 2[i ~ —0.3. 
The corresponding cumulative distribution of ./V in the tail of high values (see equa- 
tion (|19jl ) has the exponent A = 4 + 2a = 3.4. This is in an excellent agreement with 
the empirical cumulative distribution exponent 3.4 defined in [13] for 1000 stocks of 
the three major US stock markets. 

The numerical results confirm that the multiplicative stochastic model of the 
time interval between trades in the financial market is able to reproduce the main 
statistical properties of trading activity ./V and its power spectral density. The power- 
law exponents of the pdf of the interevent time, a, and the cumulative distribution 
of the trading activity A, as well as the slope of power spectral density, (3, are 
defined just by one parameter of the model 2^y/a 2 . The model suggests a simple 
mechanism of the power-law statistics of trading activity in the financial markets. 
We expect that multiplicative model of the time interval between the trades with 
more specific restrictions for the diffusion and more precisely adjusted parameters 
lies in the background of the financial markets statistics and may be useful in the 
financial time series analysis. 
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